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Abstract 

By studying the present cosmological data, particularly on CMB, SNelA 
and LSS, we find that the future fate of the universe, for simple linear models 
of the dark energy equation-of-state, can vary between the extremes of (I) a 
divergence of the scale factor in as little as 7 Gyr; (II) an infinite lifetime of the 
universe with dark energy dominant for all future time; (III) a disappearing 
dark energy where the universe asymptotes as t —> oo to a{t) pc le. 
matter domination. Our dreadful conclusion is that no amount of data from 
our past light-cone can select between these future scenarios. 


Emails: frampton@physics.unc.edu, takahash@physics.unc.edu 



Introduction 

The cosmic concordance of data from three disparate sources: Cosmic Microwave 
Background (CMB), Large Scale Structure (LSS) and High-Red-Shift Supernovae 
(SNelA) suggests that the present values of the dark energy and matter components, 
in terms of the critical density, are approximately fix — 0.7 and — 0.3. The 
question to which we try to make a small contribution in this paper is to what 
extent precision cosmological data will allow us to discriminate between possible 
futute fates of the Universe? 

If one makes the most conservative assumption that Hx corresponds to a cos¬ 
mological constant with Equation of State given by a constant w = p/p = —1, then 
the future evolution of the universe follows from the Friedmann equation 
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in which a(t) is the scale factor normalized at the present time as a(to) = 1, p(t) = 
p{0)a~^ is the energy density of matter component and A is constant. Here we 
assumed that the universe is flat and neglected radiation. 

In such a simple, and still viable, case the behavior of a{t) for asymptotically 
large t ^ oo is 


a(t) ~ exp 



( 2 ) 


so that the dark energy asymptotically dominates and the universe is blown apart 
in an inhnite time a{t) ^ oo as t ^ oo. 

Even assuming that w is constant, however, there is a wide range of possible 
w: according to [1] the allowed values are —2.68 < w < —0.78. We do not assume 
this result but will arrive at a similar allowed range; the difference is because our 
priors are slightly different (we hx the cosmological parameters Qm = 0.3, Ha = 
0.7, Hfe = 0.02, and h = 0.65 instead of allowing them to vary.) The asymptotic 
behavior in Eq.(2) is very far from established by present data. The case tc < —1 
has the property that boosting from the dark energy rest frame to an inertial frame 
with velocity satisfying (n/c)^ > —1/w leads to a negative energy density, but this 
does not violate any law of physics^. Here we shall consider some simple models for 
w, including dependence on red-shift w{z), to illustrate how far existing data are 
from answering the question of the future fate of the universe. As we will show, for 
a model in which w varies linearly with red-shift, present data are consistent with 
extremely different futures. For examples, in one case the scale factor diverges[4, 5] 
in hnite time in just another 7 Gyr, while in another case the energy density of 
dark energy decreases eventually faster than that of matter, i.e., the dark energy 
disappears and the universe reverts to being matter-dominated with a{t) t2/3, as 
t ^ oo. 


^This violates the weak energy condition [2, 3]. 

^Gravitationally-bound systems could survive longer than tr in Eq.(5) but such systems would 
be infinitely separated from one another. 
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Constant Equation of State. 


Here we discuss the future fate of the universe in the case of the constant equation 
of state. If we assume, to begin, that w is constant then, keeping only the dark energy 
term 

= H^nxa-^ (3) 

where (3 = 3(1 + w). Most authors have discussed the case with w > —1, however 
the case with tc < —1 is also possible and discussed phenomenologically in [1, 4, 6], 
and in connection with string theory in [7]. If /? < 0, corresponding io w < —1, the 
solution of Eq.(3) diverges at a hnite time t = t*. By integrating 

roo j - rt* 

/ a^/2-i _ / dt (4) 
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one hnds that the remaining time fr before time ends tr = {t* — to) is given analyti¬ 
cally by 

2 1 

~ 3Ho V^{-w - 1 ) 

In Eq.(5), putting in = 0.7 and = 9.2 Gyr one finds for w = —1.5, —2.0 
and —2.5, respectively G = 22, 11 and 7.3 Gyr. 

In such a constant w scenario which is consistent with all cosmological data, the 
divergence of the scale factor will occur in a hnite time period of 7 Gyr (or more) 
from now. 

With respect to the Solar System, this end of time occurs generally after the 
Sun has transformed into a Red Giant, and swallowed the Earth, as is expected 
approximately 5 Gyr in the future. 


Equation of State Varying Linearly with Red-Shift. 


As a more general ansatz, we consider the model for w depending linearly on 
red-shift:^ 


w{Z) = w{0) + czeic -Z) + CC9{Z - C) (6) 

where the modihcation is cut off arbitrarily at some Z = ( > 0. We assume G < 0 
and consider the two-dimensional parameter space spanned by the two variables 
tc(0) and C. 

model with w linearly depending on red-shift is also discussed in [8] but fitting the CMB 
data was not investigated. Another parametrization of w{Z) is in [9]. 
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In order to discuss constraints on our phenomenological model, we compare its 
prediction with experimental data from SNelA and CMB. We evaluate the goodness- 
of-£t parameter as a function of tc(0) and C. For SNelA, we used a dataset 
consisting of 37 SNe from [10] with the MFCS method. To calculate the CMB 
power spectrum, we used a modified version of CMBFAST [11], For the analysis of 
CMB, we used experimental data from COBE [12], BOOMERanG [13], MAXIMA 
[14] and DASI [15]. To calculate we adopt the offset log-normal approximation 
[16] and used RADPACK package [17]. We also studied the constraint from LSS 
using the 2dE data [18], and found that it does not give severe constraint on the 
parameters tc(0) and C. 

For an illustration, we take the cosmological parameters as = 0.02, IIm = 
0.3, fix = 0.7 and h = 0.65, and the initial power spectra are assumed to be scale 
invariant in all numerical calculation in this paper. 

To set the stage, let us first use only the SNelA data to constrain the parameters 
t(;(0) and C. The result is shown for (^ = 2 in Figure 1 where the 99 % C. L. allowed 
region is the region between the two dashed lines shown. We may remark three 
distinct regions: 

(I) t(;(0) < (C — 1). In this case there is divergence of the scale factor, at a finite 
future time. 

(II) (C — 1) < tc(0) < C. Here the lifetime of the universe is infinite. The dark 
energy dominates over matter, as now, at all future times. 

(III) C < tc(0). The lifetime of the universe is again infinite but after a finite time 
the dark energy will disappear relative to the dark matter and matter-domination 
will be re-established with a{t) 


When we add the constraints imposed by the CMB data, the allowed region is 
smaller as shown in Figure 2, plotted for ( = 0.5. Such a small ( still allows all 
three future possibilities (I), (II) and (HI). For somewhat larger ( only possibilities 
(I) and (II) are allowed in this particular parameterization. 


The case C = 2 is exhibited in more detail for different values of t(;(0) and C 
in Figures 3 and 4. Figure 3 shows the variation of the transition red-shift 
where deceleration changes to accelerated cosmic expansion defined by q{Ztr) = 0. 
From the figure, we can read off that Ztr becomes smaller as C becomes more 
negative for fixed ta(0); this is because the epoch where dark energy becomes the 
dominant component of the Universe becomes later. This affects the magnitude-red 
shift relation of high-Z supernovae. In Figure 4, the magnitude-red shift relation is 
shown for the SNelA data [10] along with the prediction of our phenomenological 
model for C = 2. The magnitudes are calculated, as usual, relative to the empty 
universe Milne model with = 0, Hx = 0 and = 1. One can expect that 
high-redshift SNe would appear dimmer if C were more negative. 


3 



To return to our main point, let us assume that more precise cosmological data 
will allow an approximate determination of w{Z) = f{Z) as a function of Z for 
positive Z > 0. Then to illustrate the possible future evolutions write: 

w{Z) = f{Z)e{Z) + (/(O) + aZ)e{-Z) (7) 

In this case, the future scenarios (I), (II) and (III) occur respectively for a > (/(O) + 
1)> (/(O) + 1) > a > /(O) and a < /(O). 

Present data are consistent with a simple cosmological constant f{Z) = —1 in 
Eq.(7) in which case the divergence of the scale factor occurs for a > 0, the inhnite- 
time dark energy domination for 0 > a > — 1, and disappearing dark energy for 
a < — 1 . 

Since in practice F{Z) for Z > 0 will never be determined with perfect accuracy 
the continuation of w(Z) to future Z < 0 will be undecidable from observation as 
will therefore be the ultimate fate of the Universe. 
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Figure Captions 


Figure 1. Supernovae lA constraints for C = 2 showing the three regions (I), (II) 
and (III) with different fntnres as discussed in the text. 

Figure 2. Addition of the CMB data can permit all three regions for snfhciently 
small C, only (I) and (II) for larger (. This plot is for small ( = 0.5. 

Figure 3. .^trans, dehned in the text, is plotted as a fnnction of t(;(0) for any ( > 0.7. 

Figure 4. The SNelA data and predicted curve of our phenomenological models for 
C = 2 are plotted relative to an empty nniverse Milne model with Qm = = 0 

and hlfc = 1. 


6 

































